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By analyzing CONREP data for oiler (AO) and ammunition (AE) ships 
replenishing destroyers (DD) and attack carriers (CVA) with regard to 
day and night operation and load requirement, it is shown through Chi- 
Square goodness-of- f it tests that alongside time data can^be fitted to 
an Erlang or exponential distribution. In addition, several methods 
for estimating the shape, scale, and shift parameters for gamma and 
Erlang distributions are presented. 

Employing well known but analytically little used features of the 
replenishment operation, a computer simulation model is then formulated 
and programmed based on these distributions, and thus sensitive to 
changes in load requirements. Using the simulation as an experimental 
device, an example is run to demonstrate its use in estimating the time 
to complete the operation and conducting sensitivity analyses on the 
various load requirements. 
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I . INTRODUCTION 



The primary reason the United States Navy is capable of performing 
its regularly assigned mission as well as assist in global contingencies 
effecting United States' interests, is the ability to swiftly resupply 
combatant ships almost anywhere in the world. This is possible through 
replenishment at sea operations which render ships independent of 
costly, vulnerable land based supply facilities while effectively 
increasing combatant force strength, since ships are able to remain 
on or near their assigned stations during replenishment. 

Although very important, replenishment at sea is a hazardous 
operation. During replenishment, participant ships are extremely 
close increasing their possibility of collision and vulnerability to 
enemy attack. Even though vertical replenishment (VERTREP) using 
helicopters has helped to reduce the risk in some cases, all fuel 
transfers such as Navy Special Fuel Oil (NFSO), Aviation Gasoline 
(AVGAS), and JP-5, and most transfers of dry material such as ammuni- 
tion and stores, are necessarily done by the more dangerous connected 
replenishment (CONREP) method. Recognizing the importance of CONREP, 
the U.S. Navy is using a considerable amount of its resources to 
increase the efficiency and reduce the danger of this operation. 
Increased training of personnel, redesign or replenishment ships and 
equipment, and studies utilizing analytical and computer simulation 
models are efforts that have been and are being made in this endeavor. 

Operations Research studies of the replenishment at sea operation 
have been very few. However, McCoullough formulated two models of 

a simplified underway replenishment operation. One an analytical 
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infinite cycle model using multi-stage queueing theory and the other 
a non-cyclic computer simulation. In both models the distribution of 
service time at each replenishment ship was assumed to be exponential. 
Waggoner [7*] and to some extent Patterson \j>\ also used the assumption 
of exponential service (i.e. alongside) times in replenishment at sea 
operations. Besecker \ji\ questioned the validity of this assumption 
and analyzing CONREP data, showed that the distribution of service times 
can be fitted with a gamma distribution. Categorizing the data as to 
participants and day and night operation, he contended that the data 
would also fit the more useful Erlang distribution in each case. This 
paper continues the investigation of replenishment at sea operations 
along the lines of reference \l\ . 
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II. OBJECTIVES 



The most immediate gain in CONREP efficiency can be experienced 
through proper planning of the operation. However, studies of the 
operation have been largely ignored in the planning stage, since they 
lack the realism necessary to gain the confidence of operational 
commanders. Previous models (such as \W] and U73 ) have not taken 
full advantage of the many known features of CONREP, for instance, 
arrivals at supply ships are not random but may be promulgated in 
advance by the Task Group Commander. Also, the quantity required by 
each combatant ship is known prior to the operation to some degree of 
accuracy rather than being a random variable. If the transit time 
between various supply ships is considered deterministic and known, 
the uncertainty in the system is a result only of the unknown alongside 
times, which should be dependent on the participants and the loads 
transferred. 

The first objective of this paper is to develop a more accurate 
estimation of alongside times for various combinations of CONREP 
participants by considering a given load requirement. Such conditional 
distributions would be applicable in analytical models, in computer 
simulations, and possibly in planning replenishment at sea operations, 
with more accurate predictions resulting in each case. 

The formulation and programming of a simulation model based on 
these alongside time distributions, and thus sensitive to changes in 
load requirements is the second objective. While incorporating the 
above known features, the simulation should be designed to enable its 
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user to easily adjust the arrangement of the ships involved, load 
requirements, and distribution parameters. In this regard, he could 
utilize it to test the suitability of various arrangements of combatant 
ships and the sensitivity of load changes. 



III. SOURCE OF DATA 



Data used to obtain the distribution of alongside completion times 
were from replenishment at sea (RAS) reports 3180/1A and 3180/IB for 
the period December 1967 to April 1969. Commander Service Force, U.S. 
Pacific Fleet requires these reports be made to him immediately 
following underway replenishment by those replenishment ships under 
his cognizance. Each report includes the replenishment ship type 
and hull number, receiving ship type and hull number, total alongside 
time, quantity transferred in barrels or short tones as appropriate, 
UNREP speed, whether day or night, whether EASTPAC or WESTPAC, and 
whether or not there was simultaneous port and starboard UNREP in 
progress. Full instructions concerning the submission of these 
reports are available in COMSERVPAC INST 3180. 3F of 6 November 1967. 

Alongside time is defined as the elapsed time from commencing 
approach from 300 yards astern to time of departure. Only data 
submitted by oiler (AO) and ammunition (AE) type ships involving 
destroyer (DD) and attack carrier (CVA) recipients were used, since 
a sufficient quantity was available only for these combinations. 
However, sufficient data should become available in subsequent years 
to make similar analyses possible for all combinations of participants. 
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IV. METHOD OF DATA ANALYSIS 



In order to accomplish the objective of a more accurate estimation 
of alongside time distributions for various participants, data was 
first categorized as to type of ships involved, then separated with 
regard to day or night operation, and finally subdivided into sets 
based on quantity transferred. The distinction based on day or night 
conditions was considered necessary, since from personal experience 
and as shown by Besecker, night operations consistently take longer 
than those conducted in daylight hours. The load intervals were 
selected on the basis of quantity of data available and uniformity 
of interval ranges. The latter was desirable in order to establish 
load intervals necessary for the computer simulation. 

Data in each category was analyzed by constructing histograms of 
the relative frequency with which the alongside replenishment time 
fell within ten minute intervals. The last interval is always a half 
open interval from its selected lower bound to infinity. In each case 
the histogram suggested that the distribution of the completion times 
for the respective load intervals could be fitted to a gamma or 
exponential distribution presented in Appendix A. This was expected 
since the distribution of alongside times over all intervals had been 
fitted to a gamma distribution in Ll~}. 

Figure (1) is the histogram constructed for the case in which an 
oiler (AO) was transferring 1000-1099 barrels of fuel to a destroyer 
(DD) during daylight hours. This is presented as a typical example. 
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Figure 1. Histogram of the relative frequency with which the alongside replenishment time 
fell within ten minute intervals for AO vs. DD - day, 1000-1099 barrels. 
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Alongside time (minutes) 



It was appropriate to apply a Chi-Square goodness-of-f it test to 
the null hypothesis 

H : P, = P i - (l,...,k) 
o l io 

against all alternatives, where P^ is the probability of finding a 
completion in the i^ interval with respect to the distribution in 
question, P^ o is the relative frequency of the data falling within the 
same interval, and k is the number of intervals. The probability P^ 
was calculated using Besecker's computer subroutine GAMDIS with the 
gamma parameters estimated from the sample mean and variance using 
the relationships 



-2 




where T is the sample mean and s the sample variance. 

The Chi-Square test consisted of calculating the statistic 

Q = \ 

Z ^io 

i= 1 

where is the number of datum which fell within the i^ ten minute 
interval and n is the sample size. In accordance with statistical 
theory 3 , if is true, then the random variable Q has an approxi- 
mate Chi-Square distribution. The hypothesis was accepted if 

2 

Q < X (-95, k-3) . 

The degrees of freedom was necessarily k-3 since two parameters were 
estimated from the data. Many authors caution that nP^ , i=l,2, ... ,k 
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should be greater than five, since Q has only an approximate Chi-Square 
distribution. This warning was strickly adhered to, which necessitated 
combining two or more time intervals when the quantity of data in a load 
interval was small. 

In the cases involving attack carriers, the hypothesis was accepted 
for all load intervals. Since the Erlang distribution is more ameanable 
to analytical and computer simulation models and only differs from the 

gamma distribution in that r is restricted to be a positive integer, each 

a A 

r was rounded to the nearest positive integer, the parameter “X adjusted 

according to the expression ^ 

"a - jl 

T 



and a similar test conducted. Again, the hypothesis was accepted in 
each case. Results of these latter tests are given in Tables I, II, 
III, and IV in Appendix B. Figure (2) shows a typical example with the 
probability density function superimposed on the relative frequency 
histogram for a case involving AO-CVA participants. The figure illus- 
trates the method of combining intervals consistent with the greater 
than five rule given above. 

When the above method was used for cases involving destroyer type 
ships, the hypothesis was rejected in each case. Re-examination of the 
histograms for these cases suggested that a shift of the time axis by 
an amount a would lend to an acceptance of the gamma distribution. 
Several approaches were taken to estimate the distribution parameters 
r, , and a. Although many were not successful in this study, all are 
presented for reference in future data analyses. 

Noting that the first three centeral moments of the gamma distribu- 
tion given by 
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Figure 2. Probability density function suggested by the method of moments superimposed 
on the respective relative frequency histogram for a typical AO vs. CVA case. 
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( 1 ) 



E l_T"\ = I + a, 

E (T -/c) = ^ , 

and 

E (T - O 3 = 2 I_ 

A 3 



could be used to estimate the parameters, the following relations were 
derived: 



A 

r 




* 



A 



■X = 




A 

a 



= T - 




where T is the sample mean given by 

n 



T = 1 



n 



V, T 

1=1 



i 
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( 2 ) 

(3) 

(4) 



s is the sample variance given by 

n 




and m the third centeral moment given by 

n 



m„ 



= i I 



_ 3 

(T - T) . 



i = l 



By rounding the shape parameter r to the nearest positive integer and 
adjusting the scale and shift parameters, the following expressions 
were obtained: 
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and 




where (x)* is the nearest integer to x. Using this method in several 
cases, r^ was found to be unity in each case suggesting an exponential 
distribution. Figure (3) shows a typical case with the exponential 
distribution superimposed on the histogram. The hypothesis was rejected 
for these cases as clearly suggested by the figure. 

Concern over the poor results and the fact that observed times were 
less than a, prompted another investigation into the method of obtaining 
the shift parameter. Hypothesizing an exponential distribution and 
constructing the likelihood function 



L(* f a) 



-^n e -7vnT+n> a 



V 



for a - min T 

i 



the maximum likelihood estimator for a was obtained as 



a « min T. . 
\ i 1 



This implied that 
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T- a 
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Figure 3. Suggested exponential distribution superimposed on the respective relative 
frequency histogram for a typical AO vs. DD case. 
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Alongside time (minutes) 



from Eq. (1). However, for all cases examined using this method, the 
hypothesis was rejected on the basis of the Chi-Square goodness-of- f it 
test as illustrated in Figure (4). 

Concluding that rounding to obtain an Erlang distribution was 
inappropriate when 0- ri 1, an attempt was made to use the maximum 
likelihood estimate for a under the assumption of the gamma distribution. 
The likelihood function in this case is 



Tn n 



L(a) = U ^ (T. 



rl <V> 

r- 1 j _ I A 

- a) e for ai 



min Tj 

V.* 'W 



otherwise. 



resulting in the maximum likelihood estimate of a as 



a = min T 



» 



again. Then referring to the expressions obtained using the method of 
moments and adjusting the other parameters to a, the following equations 
were derived: 

o ^ 




S 



Gamma distributions suggested by this third method did not fit the data 
as shown in Figure (5). 

In another attempt to use the method of moments to estimate the 
distribution parameters, no adjustment was made to the estimates 
given by Eqs . (2), (3), and (4). But this also proved unsatisfactory 
as illustrated in Figure (6). 

Convinced that the gamma distribution was appropriate from the 
shape of the histograms based on ten minute intervals from zero, 
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Suggested exponential probability density function superimposed on the respective 
relative frequency histogram for a typical AO vs. DD case with a estimated using 
the maximum likelihood function. 
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Figure 5. Suggested gamma probability density function superimposed on the respective 
relative frequency histogram for a typical AO vs. DD case with a estimated 
from the maximum likelihood function. 
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Alongside time (minutes) 



Figure 6. Suggested gamma probability density function superimposed on the respective 
relative frequency histogram for a typical AO vs. DD case with all parameter 
estimates unadjusted. 
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extensive probability tables were generated using Besecker's subordinate 
GAMDIS. Then by scanning these tables for similar patterns of the 
relative frequency ratios of time intervals compared to the probability 
that a completion time would be within the same interval, an Erlang 
distribution was selected that appeared to fit the data. A Chi-Square 
goodness-of- f it test was conducted and the hypothesis H accepted if 

2 

Q JC (.95,k-l) . 

This method proved quite fruitful, resulting in the acceptance of 
the hypothesis associated with an Erlang distribution in each case. A 
typical case is shown in Figure (7). Results of these tests are given 
in Tables V, VI, VII, and VIII in Appendix B with T calculated using 
the expression 




f\ 



~h 

The estimate a was chosen as the lower bound of the first interval in 
which data was observed. Figures (12) through (56) in Appendix C 
illustrate the fit of the probability density functions and the histo- 
grams for each case involving destroyer type ships. 



28 



Figure 7. Suggested gamma probability density function using scanning method imposed on the 
respective relative frequency histogram for a typical AO vs. DD case. 
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V. THE COMPUTER SIMULATION MODEL 



As stated in Naylor ^5^ the purpose of simulating a large operation 
such as replenishment at sea is to conduct situational experiments that 
would ordinarily be too expensive and/or cumbersome to perform physi- 
cally. Because of the many conceivable combinations of participant 
ships and variations in quantities transferred between them that one 
might want to examine, physical experimentation would be unreasonable 
here. Therefore, based on the distributions of the previous section, 
a Monte Carlo type digital computer simulation of the replenishment at 
sea operation named SIMA was devised and programmed in accordance with 
the second objective. 

SIMA was designed to accommodate a system of from one to four 
replenishment ships, composing an Underway Replenishment Group (URG) , 
with each ship having one or two replenishment stations on the port or 
port and starboard sides. A maximum of twenty combatant ships can be 
in the queue at each station. Furthermore, six types of combatant 
ships can be in the system. This capacity is well beyond what one 
would normally expect during an actual replenishment operation, 
therefore overload of the model is not expected. 

Drawing from personal experience, the model was formulated as 
realistically as possible, thereby taking advantage of the known 
factors of a replenishment at sea system. All combatant ships are 
Initially distributed among the various stations of the supply ships, 
with each individually characterized by the load it requires from each 
replenishment ship. Combatant ships proceed sequentially through the 
system and upon arrival at a replenishment ship it enters the smaller 
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(if more than one exists) queue. Each ship receives the services of 
every supply ship and, having done so, is ejected from the system. 
Although the combatant ships do pass through the system in a specified 
manner from their initial positions, these conditions closely resemble 
the ability of a Task Group Commander to formulate the operation prior 
to its commencement, utilizing his knowledge of the approximate loads 
transferred among the participants. 

Using the simulation as an experimental device, one is able to 

specify the participants, number of stations, initial sizes of the 

station queues, and quantities transferred between the respective 

ships, and then periodically observe the state of the system. These 

observations are made at variable time increments which are the random 

alongside completion times generated as per Appendix D. An observation 

is made and the system advances whenever any one of the URG ships has 

completed a replenishment. Each observation consists of the time since 

the start of the operation and the state of the system. The state of 

the system at these times is the number of ships, including the one 

being served, at each replenishment station and is represented by an 

n dimensional vector, where n is the number of stations in the system. 

The i*^ component of the state vector being the number of ships at the 
th 

i station. 

This model was programmed in Fortran IV and tested on the IBM 
model 360-67 digital computer at the U.S. Naval Postgraduate School, 
Monterey, California. A printout of the program is presented in 
Appendix E. Because of the digital nature of the simulation, ship 
types were given the following identification numbers: 
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Oiler - 1 

Ammunition ship - 2 
Destroyer - 1 
Attack carrier - 2 

Since all the parameters derived in the previous section are given in 
data statements, they may be revised by changing the statements. In 
using this technique, a user need only supply the following information 
to conduct replenishment at sea experiments: 



1. Type and sequence of supply ships 

2. Day or night operation 

3. Type, initial position, and load requirement of the combatant 
ships . 

This information is provided by a data deck. The first m(m - 1,2,3, or 
4) cards represent the m supply ships, each card having the following 
format: 



Columns 1-3 
Columns 4-5 
Columns 6-7 
Columns 8-10 
Column 11 



Columns 12-80 



Type of replenishment ship, i.e. AO, AE 
Number of stations available. (1 or 2) 

Supply ship identification number 
Blank 

Must be either 0,1, or blank. 

0 or blank implies that there is at least 
one other replenishment ship in the system and 
another card must be read. 

1 implies no more supply ships are in the 
system and the program continues to the next 
statement . 

Blank 



The order in which these cards appear designates the sequence in which 

combatant ships pass through the system. The m+l st card specifies day 

or night operation and has the following format: 

Column 1 Blank 

Column 2 Must be either 1 or 2. 

1 designates day 

2 designates night 

Columns 3-80 Blank 



The remaining cards give information concerning the combatant ships, 
with the cards arranged as these ships initially appear in the system. 
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The remaining cards give information concerning the combatant ships, 
with the cards arranged as these ships initially appear in the system. 
The initial condition of each station’s queue for each supply ship is 
given in accordance with the order of the supply ship cards, with the 
port station designated first. Each card represents a combatant ship 
with information given in the following format: 



Columns 1-3 
Columns 4-5 
Columns 6-11 
Columns 12-17 
Columns 18-23 
Columns 24-29 
Column 30 



Columns 31-80 



Ship type, i.e. DD, CVA 
Ship identification number. 

Load required from first supply ship. 
Load required from second supply ship. 
Load required from third supply ship. 
Load required from fourth supply ship. 
Must be either 0, 1, or blank. 

0 or blank indicates that there is 
at least one other combatant ship in 
the queue for the station in question 
and another card must be read. 

1 indicates that no more ships are 
initially in the queue in question. 
Blank 



Figure (8) shows a hypothetical initial arrangement of the supply 
and combatant ships, illustrating the possible routes of the combatant 
ships. As an example, this system was simulated using SIMA in day and 
night operation. Figure (9) shows the data cards submitted to set up the 
system for the daylight case. Since only two supply ships are in the 
example, the results may be represented graphically as in Figures (10) 
and (11). The origin represents the starting time of the operation. 

Each dot represents a combatant ship having completed replenishment at 
either the AO or AE. The length of the line segment between two 
consecutive dots indicates the elapsed time in minutes between two 
consecutive completions of the replenishment operation. For example, 
in Figure (10) a ship departed the AO at approximately 60 minutes, a 
ship left the AE at approximately 80 minutes, a ship left the AO at 
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Figure 8. Hypothetical initial arrangement of supply and 
combatant ships. 
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DD01 1460 3500 
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f DD01 2100 3010 




Figure 9. Data deck used to set up example problem. 
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Figure 10. Graphical representation of simulation 
results for the example system during 
daylight operation . 
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Figure 11. Graphical representation of simulation results for the 
example system during night operation. 
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approximately 100 minutes, a ship left the AE at approximately 100 
minutes, and so on, regarding the commencement of the replenishment 
operation 0 hours and 0 minutes. Knowing the initial arrangement of 
the ships, one is able to reconstruct the movement of each combatant 
ship through the system and its departure time from the supply ships. 

Upon running the simulation fifty times for the example, it was 
found that the expected (average) time to complete the replenishment 
operation during daylight hours was 598.62 minutes while night operations 
required 635.38 minutes. Figures (12) and (13) are the histograms for 
the frequency ratios of total replenishment time of all combatants for 
these experiments. 

We see from the data cards that the destroyer, initially the second 
ship in the queue at the port side of the AO, requires 648 barrels of 
fuel oil and 2200 short tons of ammunition. When this was changed to 
1100 barrels and 3000 short tons and the simulation rerun, the expected 
time to complete the operation was 605.54 for day and 638.10 for night 
conditions. Sensitivity analyses along these lines could easily be 
accommodated by the simulation model. 

SIMA, might also be useful in planning replenishment at sea 
operations. By tracing the path of particular ships in simulation, 
one could detect possible bottlenecks in a proposed replenishment 
system. Then by repositioning one or more ships in the formation and 
conducting additional experiments, one could plan a more efficient 
arrangement to be used in actual replenishment operations. 
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Frequency ratio Frequency ratio 



Figure 12. Relative frequency histogram for total 

replenishment time of all combatants from 
fifty simulation runs for day operation. 
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Figure 13. Relative frequency histogram for total 

replenishment time of all combatants from 
fifty simulation runs for night operation. 
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VI. SUMMARY 



In this paper we have shown that an Erlang or exponential distribu- 
tion may be fitted to the CONREP alongside time data involving AO, AE, 
DD, and CVA type ships for each of several load intervals. It is 
possible for these distributions to be used in analytical models 
employing Laplace transforms that would be sensitive to changes in 
load requirements of the combatant ships. Furthermore, several methods 
for obtaining the parameters for the gamma and Erlang distributions are 
presented for reference in future data analyses. 

The second part of the study utilizes the distributions in a compu- 
ter simulation model. Using the simulation as an experimental device, 
one can easily and realistically investigate the suitability of various 
CONREP situations confident that the results will reflect one's know- 
ledge of the approximate loads required by each of the combatant ships. 
It was demonstrated that SIMA could be used to estimate the time to 
complete the operation and investigate the sensitivity of certain load 
requirements, and suggested that the model be used in planning underway 
replenishment formations. 
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APPENDIX A 



The Gamma, Erlang, and Exponential Distributions 
The gamma probability density function is given by 



f (x) - 



T. y-v * 

n(n) 

0 



for xZ 0 
otherwise , 



where r is the shape parameter and > the scale parameter. In shifting 
the function by an amount a the random variable x would be replaced by 
x-a in the above expression. 

The Erlang and exponential distributions are special cases of the 
gamma distribution. If r is a positive integer, the probability 
density function may be written as 



f(x) = 




- r n-1 - X 
Ax e 

(n-1)! 

0 



for x- 0 
otherwise , 



and is that of the Erlang distribution. If r is unity, the distribu- 
tion is the well known exponential distribution. 
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APPENDIX B 



Tables 





Table 


I . AO vs . 


CVA - 


day 




load int. 


T 


■X 


A 

r 


Q 


T 2 (.95,k 


0-3999 


85.0 


.096 


8 


2.615 


5.99 


4000-4999 


98.7 


.244 


24 


2.716 


3.84 


5000-5999 


115.9 


.117 


14 


3.444 


5.99 


6000-6999 


127.2 


.074 


9 


1.497 


3.84 


7000-7999 


142.2 


.050 


7 


4.651 


5.99 


8000-8999 


139.3 


.156 


22 


2.411 


3.84 


9000-9999 


126.3 


.047 


6 


1.218 


5.99 


10000- 


143.3 


.060 


4 


2.584 


5.99 





Table 


II. AO vs. 


CVA - 


night 










A 






load int. 


T 


■x 


r 


Q 


"Xjt . 95,k 


0-5999 


118.2 


.068 


8 


1.007 


5.99 


6000-7999 


134.1 


.088 


12 


1.154 


5.99 


8000-9999 


139.9 


.063 


9 


0.856 


5.99 


10000-11999 


141.8 


.136 


19 


0.732 


5.99 


12000- 


194.2 


.056 


11 


0.635 


5.99 
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Table III. AE vs. CVA - day 



load int. 


T 


Tv 


r 


Q 


T 2 (.95, 


0-9999 


74.7 


.092 


7 


2.486 


3.84 


10000-19999 


105.2 


.114 


12 


2.624 


3.84 


20000-29999 


159.9 


• 111 


18 


3.055 


5.99 


30000-39999 


189.8 


.141 


27 


2.180 


3.84 


40000- 


229.4 


.079 


18 


0.067 


3.84 





Table 


IV. AE vs. 


CVA - 


night 














2 


load int. 


T 


T 


r 


Q 


X(.95,k-3) 


Or 19999 


90.1 


.084 


8 


2.573 


3.84 


20000-29999 


162.8 


.033 


6 


0.541 


3.84 


30000-39999 


168.8 


.076 


13 


1.402 


3.84 


40000-49999 


216.7 


.113 


24 


0.649 


3.84 


50000- 


258.1 


.094 


24 


0.823 


5.99 
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Table V. 


AO vs . 


DD 


- day 










A 


A 


A 


A 




r% 


load int. 


T 


T 


■X 


r 


a 


Q 


X(.95,k-1) 


0-399 


41.3 


40.0 


.100 


2 


20 


4.558 


11.10 


4Q0-499 


44.9 


50.0 


.100 


3 


20 


2.385 


9.49 


500-599 


49.8 


50.0 


.098 


2 


30 


1.973 


9.49 


600-699 


43.1 


40.0 


.100 


2 


20 


1.104 


9.49 


700-799 


47.7 


51.4 


.140 


3 


30 


3.172 


7.81 


800-899 


52.9 


53.3 


.060 


2 


20 


7.912 


12.60 


900-999 


52.8 


55.0 


.120 


3 


30 


7.143 


9.49 


1000-1099 


49.6 


55.0 


.120 


3 


30 


4.413 


14.10 


1100-1199 


53.0 


50.0 


.100 


3 


20 


5.529 


12.60 


1200-1299 


53.4 


50.0 


.100 


3 


20 


8.167 


12.60 


1300-1399 


56.4 


63.3 


.140 


5 


30 


1.389 


7.81 


1400-1499 


54.1 


55.0 


.040 


1 


30 


7.137 


12.60 


1500-1599 


62.6 


58.5 


.070 


2 


30 


4.637 


12.60 


1600-1699 


65.3 


65.7 


.140 


5 


30 


4.189 


12.60 


1700-1799 


61.3 


60.0 


.100 


3 


30 


3.306 


9.94 


1800-1899 


64.1 


60.0 


.100 


3 


30 


1.509 


11.10 


1900-1999 


67.3 


70.0 


.100 


3 


40 


2.439 


9.94 


2000-2099 


67.4 


40.0 


.100 


3 


40 


3.073 


11.10 


2100-2199 


65.2 


65.0 


.040 


1 


40 


2.306 


9.49 


2200- 


68.3 


70.0 


. 100 


3 


40 


4.925 


11.10 



Table VI. AO vs. DD - night 



A 



load int. 


T 


T 


0-599 


65.6 


67.5 


600-699 


63.3 


67.5 


700-799 


70.2 


67.5 


800-899 


78.6 


67.5 


900-999 


67.8 


67.5 


1000-1099 


61.3 


67.5 


1100-1199 


62.1 


67.5 


1200-1299 


56.1 


65.0 


1300-1399 


67.6 


65.0 


1400-1499 


71.6 


67.5 


1500-1599 


55.9 


65.0 


1600-1699 


74.4 


75.5 


1700-1799 


78.8 


75.5 


1800- 


76.7 


75.5 



A 

r 


A 

a 


Q 


X 2 (.95,k-1) 


3 


30 


5.626 


11.10 


3 


30 


1.510 


9.49 


3 


30 


0.253 


7.81 


3 


30 


1.784 


7.81 


3 


30 


0.318 


7.81 


3 


30 


1.313 


9.49 


3 


30 


2.912 


7.81 


2 


40 


1.799 


11.10 


2 


40 


8.014 


9.49 


3 


30 


1.528 


7.81 


4 


40 


3.191 


7.81 


3 


40 


0.549 


7.81 


3 


40 


1.854 


7.81 


3 


40 


4.563 


9.49 



A 

Tv 

080 

080 

080 

080 

080 

080 

080 

080 

080 

080 

160 

080 

080 

080 
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Table VII. AE vs. DD - day 



load int. 


T 


A 

T 


A 


A 

r 


A 

a 


Q 


T 2 (.95,k 


0-999 


42.8 


46.0 


.250 


9 


10 


3.584 


9.49 


1000-1999 


63.3 


60.0 


.100 


4 


20 


1.013 


9.49 


2000-2999 


77.0 


80.0 


.080 


4 


30 


3.332 


12.60 


3000-3999 


97.8 


106.6 


.060 


4 


40 


7.927 


14.10 


4000-4999 


110.5 


105.4 


.044 


2 


60 


1.524 


11.10 


5000- 


127.5 


120.0 


.080 


4 


70 


5.642 


12.60 







Table VIII. 


AE 


vs . DD 


- n 


ight 








A 


A 


A 


A 






load Int. 


T 


T 




r 


a 


Q 


T 2 (.95,k 


0-999 


45.6 


40.0 


.100 


3 


10 


5.634 


9.49 


1000-1999 


72.3 


67.5 


.080 


3 


30 


1.075 


9.49 


2000-2999 


90.9 


87.5 


.080 


3 


50 


2.913 


9.49 


3000-3999 


104.7 


100.0 


.040 


2 


50 


3.666 


9.49 


4000- 


139.7 


145.0 


.080 


3 


70 


3.149 


12.60 
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APPENDIX C 



Figures 

This appendix contains figures illustrating the fit of the 
probability density functions to the histograms for all cases 
involving destroyer type ships. 
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Figure 14. AO vs. DD - day. 0-399 barrels 
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Figure 15. AO vs. DD - day. 400-499 barrel 



f 



o 
— m 



o 




o 






o 


cn 


CSJ 


H ii 


II 


<u <K 


< <0 



CO 



I 



I 



i 

1 



y J 



r 




o 

in 



o 

m 







o 

o 

o 



49 



100 

Alongside time (minutes) 



500-599 barrels 



o 



>N 

T3 



Q 

Q 

w 

> 

O 

< 



\D 



oj 

u 

s. 

*r4 



j 

I 




/ i 

« 

/ ! j 





00 






ON 






o 


O 


CM 


• 


CO 


II 


H 


n 


< u 


<r< 


<« 




o 

<r 



o o 

co cm 




Aouanb^a^ 



50 



Alongside tire (minutes) 



Figure 17. AO vs. DD - day. 600-699 barrels 
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Figure 19. AO vs DD - day. 800-899 barrels 
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Figure 20. AO vs. DD - day. 900-999 barrels 
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Figure 21. AO vs. DD - day. 1000-1099 barrels 
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Figure 22. AO vs. DD - day. 1100-1199 barrels 
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Figure 23. OA vs. DD - day. 1200 - 1299 barrels 
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Figure 24. AO vs. DD - day. 1300-1399 barrels 
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Figure 25. AO vs. DD - day. 1400 - 1499 barrels 
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Figure 26. AO vs. DD - day. 1500 - 1599 barrels 
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Figure 28. AO vs. DD - day. 1700 - 1799 barrels 
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Figure 29. AO vs. DD - day. 1800 - 1899 barrels 
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Figure 30. AO vs. DD - day. 1900 - 1999 barrels 
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Figure 31. AO vs. DD - day. 2000 - 2099 barrels 
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Figure 32. AO vs. DD - day. 2100 - 2199 barrels 
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Figure 33. AO vs. DD - day. 2200 - barrels 



\ 




-1 _ _ 

o o o 

CO CM r-H 



o’ptfj Aouanb 



o 



o 

o 



67 



Alongside time (minutes) 



Figure 34. AO vs. DD - night. 0 - 600 barrels 
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Figure 35. AO vs. DD - night. 600 - 699 barrels 
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Figure 36. AO vs. DD - night. 700 - 799 barrels 



I 



I 



I 



o 

m 




-l 1 _ .. JL 

O o o 

m cn h 






o 



o 

o 



70 



Alongside time (minutes) 



899 barrels 



I 



I 




I 

o 

o 

00 



si 

00 

tA 

c 

I 

o 

o 



CO 

> 

o 

< 



co 

GJ 

U 

3 

00 

•H 




o o o 

CO (N *— < 

• • • 

oji UJ X3u»nbaj,£ 



o 

o 



71 



Alongside time (minutes) 



Figure 38. AO vs. DD - night. 900 - 999 barrels 
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Figure 39. AO vs. DD - night. 1000 - 1099 barrels 
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igure 40. AO vs. DD - night. 1100 - 1199 barrels 
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Figure 41. AO vs. DD - night. 1200 - 1299 barrels 
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Figure 42. AO vs. DD - night. 1300 - 1399 
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night. 1400 - 1499 barrels 
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Figure 44. AO vs. DD - night. 1500 - 1599 barrels 
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Figure 45. AO vs. DD - night. 1600 - 1699 barrels 
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Figure 46. AO vs. DD - night. 1700 - 1799 barrels 
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Figure 47. AO vs. DD - night. 1800 - barrels 
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Figure 48. AE vs. DD - day. 0 - 1000 short tons 
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Figure 49. AE vs. DD - day. 1000 - 1999 short tons 
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Figure 50. AE - vs . DD - day. 2000 - 2999 short tons 



! 

! 




x 

O 




ojivjl Aouanbej^ 



84 



Alongside time (minutes) 



Figure 51. AE vs. DD - day. 3000 - 3999 short tons 
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Figure 52. AE vs. DD - day. 4000 - 4999 short tons 
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Figure 53. AE vs. DD - day. 5000 - short tons 
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Figure 54. AE vs. DD - night. 0 - 1000 short tons 
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Figure 55. AE vs. DD - night. 1000 - 1999 short tons 
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Alongside time (minutes) 



Figure 56. AE vs. DD - night. 2000 - 2999 short tons 
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Alongside time (minutes) 



Figure 57. AE vs. DD - night. 3000 - 3999 short tons 
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Alongside time (minutes) 



Figure 58. AE vs . DD - night. 4000 - short tons 
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APPENDIX D 



Generating Random Observations from an Erland Distribution 

The most important operation in a Monte Carlo type computer 
simulation is the generation of appropriate random variates. In this 
study, it was necessary to obtain random variates from the Erland 
distribution. Methematically the Erland distribution is a convolution 
of r independent identical exponential distributions, i.e., the 
distribution of the sum of r exponentially distributed, independent 
random variables whose density is given by 

f (y) = e' >y , Y>0. 

Thus, we are able to generate Erlang random variables by summing r 
independent exponentially distributed variables , Y , ... , Y r , 

each with expected value 1/ . Doing so we obtain that 

y i 

i=l 



has the Erland distribution with parameters r and , i.e. with density 



f(x) 



"X 



r r-1 ->x 

x e 

(r-1) : 

0 



for x 1 0 
otherwise . 



The question now is how to obtain exponentially distributed random 
variates? The cumulative distribution function for the exponential 
distribution is given by 
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f(y) = l-e'^ 7 , y-0. 

It is well known that F(Y), and by symmetry l-F(Y) also, is a uniform 
random variable over the interval (0,1) regardless of what the distri- 
bution of the continuous random variable Y is. This is equivalent to 
setting 

-AY 

W = e 

where W is a random variable from the uniform (0,1) distribution. 
Therefore, taking the natural logarithm of both sides 

-AY = In (W) 

which yields 

Y = - — ln(W) . 

* 

Thus, if W is a uniform random variable over (0,1), Y will have the 
desired exponential distribution. Substituting into Eq # (1) we get 

X =- -Lln(TT W i ). 

A i=1 

There are several methods available for generating uniform (0,1) 
random variables. The IBM Scientific Subrouting RANDU was used in 
the study. 
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APPENDIX E 



Computer Program 



This appendix contains a printout of the computer program SIMA. 
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